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We show how a quantum Ising spin chain in a time-dependent transverse magnetic field can be 
simulated and experimentally probed in the framework of circuit QED with current technology. The 
proposed setup provides a new platform for observing the nonequilibrium dynamics of interacting 
many-body systems. We calculate its spectrum to offer a guideline for its initial experimental 
characterization. We demonstrate that quench dynamics and the propagation of localized excitations 
can be observed with the proposed setup and discuss further possible applications and modifications 
of this circuit QED quantum simulator. 
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The promising idea of tackling complex quantum 
many-body problems by quantum simulations [TJ |2] 
has become even more compelling recently, due to the 
widespread current interest in nonequilibrium dynam- 
ics. Indeed, experiments with cold atoms in optical lat- 
tices [3-6 J and ions [7HT0] have already made impressive 
progress in this regard. At the same time, the capabilities 
of scalable, flexible solid-state platforms are developing 
rapidly. In particular, circuit quantum electrodynamics 
(cQED) architectures of superconducting artificial atoms 
and microwave resonators [TTHT5] are now moving to- 
wards multi-atom, multi-resonator setups with drasti- 
cally enhanced coherence times, making them increas- 
ingly attractive candidates for quantum simulations [16] . 
Here we propose and analyze a cQED design that sim- 
ulates a quantum transverse-field Ising chain with cur- 
rent technology. Our setup can be used to study quench 
dynamics, the propagation of localized excitations, and 
other nonequilibrium features, in a field theory exhibit- 
ing a quantum phase transition (QPT) [17], and based 
on a design that could easily be extended to break the 
integrability of the system. 

The present work takes a different path than the 
proposals for simulating Bose-Hubbard-type many-body 
physics in arrays of coupled cavities, which might be also 
realizable in the framework of cQED p~6| fT8| [T9] , It is 
based on an independent, possibly simpler concept - the 
direct coupling of artificial atoms - that naturally offers 
access to quantum magnetism. The transverse-field Ising 
chain (TFIC) is a paradigmatic quantum many-body sys- 
tem. It is exactly solvable j20[|2T] and therefore serves as 
standard theoretical example in the context of nonequi- 
librium thermodynamics and quantum critical phenom- 
ena [T7l [22---27J . Our proposal to simulate the TFIC and 
its nonequilibrium dynamics might help to mitigate the 
lack of experimental systems for testing these results. 
Moreover, the experimental confirmation of our predic- 
tions for various nonequilibrium scenarios in this inte- 
grable many-body system would serve as an important 



benchmark and allow the experimentalist to proceed to- 
wards variations of the design that break integrability or 
introduce other features. 

Implementation of the TFIC. — A charge-based ar- 
tificial atom (such as the transmon) [28J in a supercon- 
ducting microwave resonator can be understood as an 
electric dipole (with dipole operator a x ) that couples to 
the quantized electromagnetic field in the resonator . 
Consider the system shown in Fig. [I] (at first, without res- 
onator B). Only the first artificial atom couples to res- 
onator A. However, all atoms will couple directly (not 
mediated by a quantized field) to their neighbours via 
electric dipole-dipole coupling oc o- l x a 3 x . Since this inter- 
action is short-ranged, we model the system by 

U = uj (a^a + 1/2) + g(o) + a)al + Hi, (1) 

where T~Lj is the Hamiltonian of the TFIC, 

Q N N ~! 

Ki = ^J2c-i-jj244 +1 . (2) 

j=l 3 = 1 

Here, a) generates a photon with frequency c^o, and (J 3 x j z 
is a Pauli matrix. That is, we consider the artificial atoms 
as two- level systems (qubits). This is justified even for 
weakly anharmonic transmons since the experiments pro- 
posed below involve only low atomic excitation probabil- 
ities or well controllable excitation techniques (7r-pulses) . 
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FIG. 1. Circuit QED implementation of the Ising model with 
transverse magnetic field. The dipole moments of the artificial 
atoms tend to align. Resonator A (B) facilitates initialization 
and read-out of the first (iVth) artificial atom by standard 
circuit QED techniques. 
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Qubit 1 and the resonator couple with strength g. The 
qubit level spacing Q > is tunable rapidly (< Ins) 
[HJ H21 [TH [15] and corresponds to the transverse mag- 
netic field in the usual TFIC. In our geometry, the qubit- 
qubit coupling strength J is positive (ferromagnetic; anti- 
ferromagnetic coupling J < can be achieved by rotat- 
ing each qubit in Fig. [I] by 90° and is discussed in |33j). 
Estimates based on the typical dimensions of a cQED 
system yield J/2tt ~ 100MHz. Interdigitated finger ca- 
pacitors between the qubits might significantly increase 
J. Resonator A facilitates initialization and read-out the 
first qubit (with standard cQED techniques, [H]). Res- 
onator B would allow one to measure end-to-end correla- 
tors. However, for simplicity we consider a system with 
only one resonator unless noted otherwise. We mention 
that the proposed setup should be also implementable us- 
ing the novel, high-coherence 3d cQED devices [30J. For 
a different proposal on the implementation of and mean- 
field-type experiments with the TFIC in cQED, see |31| . 

In our calculations (|33J) we frequently use the spin- 
free- fermion mapping for Hi from Refs. [20) 121 ] . It yields 
Hi = ^Zf, ^k(VkVk — 1/2), where rj\ generates a fermion 
of energy = 2J^/ 1 + £ 2 — 2£ cos and £ = Q/2 J is 
the normalized transverse field. The possible values of k 
are solutions of sinkN = £ > sink(N + 1). For N — >• oo, 
Hi undergoes the second order QPT at £ = 1 from a 
ferromagnetic phase (£ < 1) with long-range order in 
<j x to a disordered, paramagnetic phase (for details, see 
[IB [201 [22 [33]). 

Spectrum of the system. — An initial experiment 
would likely seek to characterize the setup by measuring 
the transmission spectrum S of the resonator as a func- 
tion of probe frequency uj and transverse field £. To cal- 
culate 5, we first focus on the spectrum of the bare TFIC, 
p(uj) = J dte luJt (o~l(t)al(0)) . It shows at which frequen- 
cies a field coupled to o\ can excite the chain. Assuming 
g/uoo <C 1, we may then approximate the chain as a lin- 
ear bath, coupled to the resonator. That is, we replace it 
by a set of harmonic oscillators with the spectrum p(uj) 
of the TFIC, which allows us to compute the resonator 
spectrum S. We calculate all expectation values for zero 
temperature. Except for near-critical coupling where Hi 
becomes gapless, this is experimentally well justified. 

For finite N, the calculated spectrum p{uj) would con- 
sist of discrete peaks. In an experiment, they would be 
broadened by decay and, for sufficiently large TV, the 
measured spectrum would be continuous. This can be 
modeled by taking N — >• oo in our calculations. In that 



p(u) = 2nS(u)G(l-0(l-e) 
+ ^R ev /1 - cos 2 jfc(w) 



(3) 



for uj > 0, and p(uj < 0) = 0. Here, Q(x) is the Heaviside 
step function, and cosk(uu) = [1 + £ 2 - (^j) 2 ]/2£. The 



delta function for £ < 1 is due to the finite mean value 
of Re(al(t)al(0)) in the ferromagnetic phase. We plot 
p(uj) for several £ in Fig.|2|a). For £ > 1 (£ < 1), p has a 
width of 4 J (4J£), which agrees with the bandwidth of 
Afc. This might be helpful to measure J. At the critical 
point £ = 1, p becomes gapless and, thus, carries a clear 
signature of the QPT. The loss of normalization for £ = 
0.5 is compensated by the delta function in (J3|. This 
is required by a sum rule for p and can be understood: 
In the ordered phase, the ground state |0) of the TFIC 
becomes similar to a a^-eigenstate. Thus, driving via a\ 
is less efficient in causing excitations out of |0), but a 
static force on o\ will change the energy of |0). We note 
that for all £, p(u) has its maximum where the band A& 
has zero curvature (and maximum slope). Thus, most rj^- 
excitations of the TFIC have a nearly uniform velocity 
= max[dAfc/d/c] (vo = 2J£ for £ < 1 and vo — 2 J for 
£ > 1), which will be important below. 

With this result for p((j), we obtain the spectrum S(u) 
of the resonator, 



S(u>) = 



[w 2 /w - w - Vx(w 2 )] 2 + [k + g 2 p(oj)} 2 ' 



(4) 



Here, x( Co ' 2 ) denotes the principal value integral x(^ 2 ) = 
1/(2tt) Jdftp(n)n/(uj 2 - n 2 ) and k is the full linewidth 
at half maximum of the Lorentzian spectrum of the un- 
coupled (g = 0) resonator. Our calculation makes use of 
tools explained, e.g., in [32], and is actually general, i.e., 
it also applies when the resonator couples to a different 
system, with another spectrum p(oj). We plot S as func- 
tion of uj and £ in Fig. |2|b). For comparison, we also 
plot the resonances of the well-known Jaynes-Cummings 
model, as they have been observed in numerous cQED 
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FIG. 2. Spectrum of the system, (a) Spectrum p(uj) = 
f dte luJt {al(t)al(0)) of an isolated transverse- field Ising chain 
for N — »> oo and normalized transverse fields £ = Q/2J — 
4, 1.2, 1, 0.5. (b) Spectrum S of a resonator coupled to a TFIC 
(as in Fig.[T]), plotted vs. £ and uj (for N —> oo). The param- 
eters used are g = 0.12, J = 0.1, and k = 10 -4 (in units 
of ujo). For better visibility of the features, values > 3.8 are 
plotted white. Dashed, excitation energies of H for N — 1. 
(c) S vs. uj for £ = 3.9,4.85,6.1 (blue, red, green). These 
lines correspond to cuts along the arrows in (b). 
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experiments (dashed; special case N — 1 in T~L). As long 
as the spectrum p(uS) of the chain does not overlap the 
resonator frequency uoq , there is a dispersive shift in close 
analogy to the off-resonant single-qubit case. Here, the 
chain is visible in S(uo) only as a small but broad side 
maximum and hardly modifies the dominant Lorentzian 
peak [green and blue lines in Fig.[2jc)]. If the chain comes 
into resonance, this changes dramatically, and S(u) takes 
on large values over a region of width ~ 4J. For our 
choice of parameters, S(u) develops a slightly asymmet- 
ric double-peak structure [red line in Fig. [2jc)]. This is 
again reminiscent of the Jaynes-Cummings doublet, but 
now the peaks are split by 4 J rather than 2g. We empha- 
size that the shape of the spectrum on resonance depends 
significantly on the ratio g/J. The larger g/J > 1, the 
closer the system resembles the single-qubit case, which 
corresponds toJ = 0. If gj J < 1, the double-peak struc- 
ture vanishes and one observes a Lorentzian around uq 
with full width at half maximum 2g 2 / J (for g 2 / J ^> k). 
This is because the resonator irreversibly decays into the 
chain, whose inverse bandwidth oc 1/J sets the density 
of states at uo « ujo and thereby determines the decay 
rate (for plots on both limiting cases and for finite TV, 
see [33J). 

Propagation of a localized excitation. — If the trans- 
verse field £ is chosen such that chain and resonator are 
off-resonant, they are essentially decoupled. In this situ- 
ation, the proposed setup allows one to study nonequilib- 
rium dynamics in the TFIC. The resonator can be used 
to dispersively read out the first qubit. For measure- 
ments, this qubit must be detuned temporarily from the 
rest of the chain such that it dominates the dispersive 
shift of the resonator [11 J and becomes decoupled from 
the chain's dynamics. 

First, we focus on the nonequilibrium dynamics of the 
chain after a local excitation has been created. As the 
resonator couples only to one qubit, the initialization of 
the system is easy. We assume that the chain is far in 
the paramagnetic phase (£ ^> 1). Hence, (a J z ) « — 1 in 
its ground state. By applying a fast (~ Ins) 7r-pulse, the 
first spin of the chain can be flipped without affecting the 
state of the other qubits (if J/2tt <C 1GHz/27t, or if the 
first qubit is detuned from the others for initialization). 
We model the state of the system immediately after the 
7r-pulse by <r^|0), where |0) is the ground state of the 
TFIC. The time evolution of the qubit excitations (cr J z ), 

(ai)(t) = (Olale^aie-^WJO), (5) 

is plotted in Fig. [3] for a chain with N = 20 and £ = 8 
(right panel). The experimentally measurable trace of 
(o-\)(t) is singled out left. Due to the qubit-qubit cou- 
pling, the excitation propagates through the chain, is re- 
flected at its end, and leads to a distinct revival of (<r\) at 
Jt R w N. Assuming J/2tt = 50 MHz, we find t R « 64 ns 
for N = 20, which is safely below transmon coherence 
times. Note that the excitation propagates with velocity 



vo = 2J. This is because a localized excitation in position 
space is a superposition of many excitations in /c-space, 
and we have seen that most of them have velocity vq. 

Quench dynamics. — One of the most appealing ap- 
plications of the proposed system would be to observe its 
nonequilibrium dynamics after a sudden change of the 
qubit frequencies Vt and, thus, of the transverse field 
£ = Q/2J. By using fast flux lines, such a change 
can be achieved practically instantaneously on the dy- 
namical time-scale (without changing the wave function) 
H3H1HS]. This (global) quantum quench in the TFIC 
has been studied theoretically, e.g. in [23-26J. One usu- 
ally assumes that for t < the system is in the ground 
state \0) a of the Hamiltonian Hi ia a t some initial value 
£ a = Q a /2J. At t = 0, f is changed to £ b = Q b /2J ", 
and the time evolution under the action of is inves- 
tigated. 

In the following, we focus on the dynamics of the exper- 
imentally easily accessible observable after quenches 
within the paramagnetic phase. This corresponds to our 
estimates for realistic values of J. The main difference of 
quenches involving the ferromagnetic phase would be the 
different dynamical time scale due to the the different 
value of vq in this phase. Fig. [4] shows the magnetiza- 
tion (cr J z )(t) after quenching £ in a TFIC (center). In re- 
gion I (see schematic plot, right), the magnetization first 
increases and then oscillates with decreasing amplitude. 
Here, it is virtually identical with the overall magnetiza- 
tion of a cyclic TFIC with N — >• oo calculated in [23J , and 
would, for N — > oo, approach a constant value. This is 
in line with predictions from conformal field theory |25| . 
However, at t = j/vo and t = (N — j)/vo (dashed red 
lines in the schematic plot), where vq = 2 J as before, the 
magnetization has dips. They are followed (in regions II 
and III) by a relaxation similar as in region I to the same 
asymptotic value (see [33J for a zoomed-in plot). Near 
the boundaries of the system, the magnetization reaches 
and stays at this value for a considerable time before un- 
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FIG. 3. Propagation of a localized excitation. Right, 
nonequilibrium time evolution of (a J z ) for all qubits j of a 
transverse Ising chain of length N — 20 in a normalized trans- 
verse field £ = Q/2J = 8 (paramagnetic phase) after the first 
spin has been nipped. Values > — 0.5 are plotted white. Left, 
separate plot of (al) on the same time scale. This quantity 
can be measured in the setup of Fig. [l] 
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dergoing a revival. A sharp oscillation across the entire 
chain at T = N/vq subsequently decays. Revivals re- 
occur (quasi-)periodically with period T (region IV), but 
this behaviour is smeared out for large times (not plot- 
ted). These phenomena are reflected in the measurable 
observable (cj\){t) (left panel) and take place on a time 
scale of - 0.1 fis for TV = 30 and J/2tt = 50 MHz. 

Our results can be qualitatively understood in a sim- 
plifying quasiparticle (QP) picture that has already been 
used to calculate or interpret the (quench) dynamics of 
different quantities in the TFIC [24-27J. In the param- 
agnetic phase, the QPs correspond to spins pointing in 
the -fz-direction. They are created in pairs by the quench 
and ballistically move with velocities with reflections 
at the boundaries. Further, only contiguously generated 
quasiparticles are correlated. After an initial transient, 
any given site will be visited only by QPs originating from 
distant places which are completely uncorrelated. This 
leads to the relaxation of the magnetization towards a 
steady-state value in region I that would be characterized 
by a certain static density of uncorrelated QPs. However, 
once correlated QPs meet again due to reflections at the 
boundaries, coherences are re-created and will show up 
in oscillation revivals. This happens, first, at multiples 
of T (black solid lines in the schematic plot) when all 
QP trajectories cross their momentum-inverted counter- 
parts (the solid red lines show an example), and, second, 
along the trajectories of QP pairs generated at the edges 
of the chain. Such QPs travel together as one partner 
is reflected at t « (dashed red lines; not plotted in 
region IV for clarity). The periodicity of the QP trajec- 
tories should lead to periodic revivals for t > T. This is 
indeed observed approximately, although finally the ve- 
locity dispersion of the QPs renders the time evolution 
quasiperiodic. Finally, QP-trajectories cannot intersect 
at j = 1, TV. The density of (incoherent) QPs is therefore 
lower here than for bulk sites, resulting in an appreciably 
lower quasistationary value. 



Discussion and outlook. — The setup and the experi- 
ments we have proposed might help to establish the sim- 
ulation of interacting quantum many-body systems as a 
new paradigm in circuit QED and to bring parts of the 
theoretical discourse in nonequilibrium physics closer to 
observation. All results discussed here should be observ- 
able within the system's coherence time, and we have nu- 
merically verified that they are robust against disorder up 
to a few percent in Q and J. Detuning individual qubits, 
however, would allow one to create arbitrary potentials 
for the excitations, study the interplay of Anderson lo- 
calization and many-body physics, or change the effective 
chain length. Using a second resonator, the dynamics of 
the end-to-end correlator (crla^) (indicating long-range 
order) could be measured (see |33j). Many other experi- 
ments are conceivable with our setup, like suddenly cou- 
pling two isolated chains (and other local quenches) or 
even parameter ramps through the QPT, with Kibble- 
Zurek defect creation. We note also that hitherto unex- 
plored measurement physics could be studied when the 
first qubit is not detuned from the chain, like resolving 
many-body eigenstates or the Quantum Zeno effect in a 
many-body system. Once the setup is properly under- 
stood, it will be easy to break the integrability of our 
model in a controlled way (e.g., via longer-range cou- 
plings). This would push our cQED quantum simulator 
into a regime beyond classical computational capabili- 
ties, where further open questions about nonequilibrium 
dynamics can be addressed, like thermalization and diffu- 
sive transport. Furthermore, going to 2d or 3d introduces 
new design options, for instance, frustrated lattices. 
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UC Berkeley for their hospitality. Support by NIM, the 
Emmy-Noether program, and the SFB 631 of the DFG 
is gratefully acknowledged. 
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FIG. 4. Behaviour after a quench: time evolution of the 
magnetization (a J z ) in a TFIC of length N — 30 after a quench 
of the normalized transverse field £ = 8 —> 1.2 (center) with 
a schematic plot (right) and the measurable observable (al) 
singled out (left) on the same time scale. Values < —0.9 
(> —0.6) are plotted black (white). 
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DIAGONALIZATION AND SPECTRUM OF THE 
TRANSVERSE-FIELD ISING CHAIN 

In this section, we diagonalize the Hamiltonian Hi 
[Eq. Q of the main text] and calculate the qubit au- 
tocorrelator p(t) = (^(^^(O)) and the corresponding 
spectrum p{uj) = j dte zu}t p(t) . Our method and notation 
follow Ref. [1J. 

In the main text, we have focussed on a circuit QED 
system with ferromagnetic qubit-qubit coupling J > 0. 
Since setups with anti-ferromagnetic coupling are also 
conceivable, we generalize in these supplementary notes 
the Hamiltonian of the transverse-field Ising chain to 



n 



N 



N-l 



+i 



(S.l) 



where J may be negative (Vt > as before). We define 
J = Applying the Jordan- Wigner transformation 



cr+ = cj exp(z7r J2k=i c l c k) to leads to 



N 



N-l 



= - ^ + « 2 c\cj -jJ2 + 4 c i+i + Hx -] ■ 

i=i j=i 

(S.2) 

with fermionic Cj. In this form, T~Li can be diagonalized 
using the method for diagonalizing quadratic fermionic 
Hamiltonians of the form 

N 

H=J2 >•!•>-,/':< + V2(4i^4 + H.c.)] (S.3) 

*.J=1 



of Ref. [lj. In our case, 



A 



( n 

-J 





V o 



-J 

-J 





-J 



-J 

-j n 
-J 



-J 



(S.4) 



and B is obtained by substituting = Q — and 
Ai+i^i = — J —> J in A. H is diagonalized by introduc- 
ing new fermions rjk = YljLi 9k,j c j + hkjCj. The com- 
ponents gkj and hkj of the vectors and hk and the 
eigenvalues A& of H are determined by defining normal- 
ized vectors <pk = 9k + hk and ipk = 9k — hk and solving 
the equations 

<t> k {A -B)= A^, M A + B) = A k (j>k. (S.5) 

For Afc 7^ 0, this is most easily done by solving, e.g., 

(A-B)(A + B)<t> k = A 2 k ct> k (S.6) 



and calculating ipk via Eqs. ( |S.5[ ). Note that since A T = 
A and B T = —B, > and the <pk and ^ can be 
chosen real and orthogonal for different fc, ^ <f>kj<l>k',j = 
J2j ^kjipk'j — Sk,k f - For A and B as defined above, one 
obtains 



W/ = ^A fc (^77 fc -l/2), 



<t>k,j 

^kj 



2JVl + £ 2 -2f cosfc, 
A k smk(N + l-j), 



sign 
2(27V 



A k sin kj, 



sin/c(7V + 1) 

l-sin[fc(2iV + l)]/sinfe) 



-1/2 



(S.7) 

(S.8) 
(S.9) 

(S.10) 
(S.ll) 



Here, f = VL/2J is the normalized transverse field, and 
the possible values of k are solutions of 



sin kN 
sinfc(iV + 1) 



(S.12) 



If |f | > N/(N + 1) (|f | < N/(N + 1)), Eq. (Sl2) has TV 
(N - 1) real solutions e [0,tt]. If |f | < A/"/ (N + 1), there 
is also one imaginary solution fc' = in (k' = n + m) for 
positive (negative) f with sinhftA/y sinh^(A/" + 1) = |f |. 
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of Eq. (S.12), 



1.0 1.5 2.0 

wave vector fc 

FIG. SI. Excitation energies A& vs. the allowed real wave 
vectors fc of a transverse-field Ising chain with N = 30 and £ = 
0.5, 1, 1.2,3 (dots). For £ = 0.5, there is also one imaginary 
wave vector (see text). The solid lines are a guide to the 
eye. The shaded regions indicate the bandwidth of the Ising 
chain for £ = 0.5 (ferromagnetic phase, orange) and for £ = 3 
(paramagnetic phase, blue). 

These solutions exhaust the eigenmodes of the system. 
Note that A fc / -> if |f | or TV -> oo. 

For A/" — >• oo, undergoes a second order QPT at 
|£| = 1 from an (anti-)ferromagnetic ordered phase (|£| < 
1) with doubly degenerate eigenstates (A^/ — > 0) to a 
paramagnetic disordered phase with A& > for all fc. It is 
signaled by correlators of the order parameter a x . Note, 
though, that (a J x ) = for all £. Since commutes with 
Y\j & J Z , all eigenstates of Hj formally obey this symmetry 
that maps a J x — > — cr^. 

Fig. [ST] shows the excitation energies A& of Hi vs. the 
allowed (real) wave vectors fc for N = 30 and various 
£ (for £ = 0.5, there is one imaginary wave vector k' « 
0.693i, and A&/ ~ 0). In the limit N —> oo, the A& form a 
continuous band. Its gap is given by |1 — 1£|| and vanishes 
at the quantum critical point |£| = 1. In the disordered 
phase (|£| > 1), the bandwidth is 4 J (indicated for £ = 3 
in Fig. SI) and independent of £. In the ordered phase 
(|£| < 1), the bandwidth is given by 4J|£|. For more 
details on the transverse- field Ising chain and its QPT, 
see, e.g., [TH3]. 

Assuming zero temperature, the qubit autocorrelator 
p(t) = (o~l(t) crl(0)} can now be easily calculated using 



c\ +ci 



k 



One obtains 



>L e ~ 



(S.13) 



(S.14) 



h>i = — arctan 

7T 



f sin(7rZ/JV) 



€cos(nl/N) - 1 



■0(1/N). (S.15) 



With these expressions for ki and ^, p(t) can be trans- 
formed into an integral j ± dl for N — >> oo. Substituting 
dZ — » dfc (fc as defined above, dfc/dZ ~ 7r/7V) and dropping 
all terms 0(1/N) finally leads to 



p(t) = 6(l-|£|)(l-|£| 2 ) 



£ 2 sin 2 fc 



dfc- 



l + £ 2 -2£cosfc 



-itA(fc) 



(S.16) 



where 6(x) is the Heaviside step function and A(fc) 
stands for A& with continuous fc. The first term on the 
RHS of Eq. pl6|> is the fc'-term in Eq. (|STl4|> for N -> oo, 



which must be treated separately. It causes a nonzero 
mean value of Rep(t) in the ordered phase. Taking the 
Fourier transform of Eq. ( S.16| ) yields 



p(uj) = 27T(5(cj)9(1 - 
+ 6(cj-2J|1 

4iei 



i^i)(i-i^i 2 ) 
-ieii)e(2j|i- 



x -^Vl-cos^H, 



(S.17) 



where cosfc(cj) = [1 + £ 2 - (|j) 2 ]/(2£). Note that this 
result does not depend on the sign of J (and £ = Q./2J). 
For ferromagnetic coupling > (and £ > 0), Eq. ( |S.17 ) 
can be simplified to the form of Eq. ([3| and is plotted 
m Fig. |2fa) of the main text. For anti-ferromagnetic 
coupling J < (and £ < 0), one just has to replace 
£ — >• |£| in Eq. (J3|. Thus, with this replacement, our 
discussion of p(cj) below Eq. Q and the plots in Fig. [5] 
of the main text hold for anti-ferromagnetic coupling as 
well. 



SPECTRUM OF THE RESONATOR 

In this section, we calculate the spectrum S(uj) of the 
resonator of our system (which is coupled to the Ising 
chain) and, complementary to Figs. [2^b) and|2|c) of the 
main text, plot S(u) for the limiting cases g/J <C 1 and 
g/J^> 1, and for finite N. 

In order to calculate S(lj), we assume g/uo <C 1 and 
linearize the Hamiltonian T~L [Eq. of the main text]. 
That is, we now consider the Hamiltonian 



The Fourier transform p(uj) of p(t) is a sum of delta 
peaks. In order to obtain a continuous spectrum p(cj), 
we have to take the limit N — >> oo in Eq. (S.14). As 
its RHS contains rapidly oscillating terms for N —> oo 
(like sin iVfc), it cannot be straightforwardly transformed 
into an integral via a Riemann sum. We therefore write 
fc/ = tt/N(1 — v\) for Z = 1, . . . , N [I] and find, by means 



N 



■ ljqXq) + x £ ^j X 3 + (S.18) 



where A. 



^2g 2 uj Xj and ?4 = SjLi(P 



: x 2 )/2. It 

is obtained by substituting the coupling term in Eq. M 
by g(at -fa) XljLi ^j x j ftj 1S a coupling constant) and 
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Hi by H^, the Hamiltonian of a set of harmonic oscil- 
lators with frequencies Wj, and by introducing canonical 
coordinates for the resonator via x$ = 1/ y/2ujo(a^ + a) 
and po = i^wo/2(cfi — a). Note that couples to a force 



Fh(f) = EjLi in E( i- P^ - B y writin g ft t E( i- 

Q ] in terms of xq and po? one finds that here xo couples 
to a force -Fj(t) = ^2g 2 ujQcrl,(t). The parameters Aj and 



Wj in % can be chosen such that 



N 



<W,(0)> = ^4- 

.7 = 1 J 



A^ . , 

J g — lWjt 



(Fj(i)Fj(O)) (S.19) 



(in this case also the spectra of the forces will agree). 
Indeed, Wj = A^. and X 2 = 2wjA 2 . sin 2 Nkj guarantee 
Eq. |Sl9| ). We now calculate 5(cj), the Fourier transform 
of 2cJo(0|xo(t)xo|0), where |0) is the ground state of T~L. 
To that end, we first reformulate 



Thus, we have expressed the spectrum of the resonator 
S(u) in terms of the spectrum (F^F^)^ of the bath of 
harmonic oscillators which is the Fourier transform of 
(F h (t)F h (0)) [Eq. flS.19D]. Note that in the limit N 



oo, (FhFh)oj can become continuous and then x(clt) is 
a principal value integral. If we now assume that we 
have chosen Aj and wj in H [Eq. ( |S. 18 )] such that Eq. 
( |S.19| ) holds, we can substitute (FhFh)^ {F I F I ) 0J = 
2g' z ujQp(uj). This leads to 



5(w) 



4e(w)w [5 2 w p(w) + 0+] 



[W 2 - W 2 -4 5 2 W0X ( W 2)]2 + [0+ + g2 LJo p(co)}l 



where x( w2 ) is the principal value integral 

. 2 , 1 I - 
X(- 2 ) = ^ / ^^2- 



(S.29) 



(S.30) 



ft 



with X T = (xq, x\, . . . , xn), P 2 

/ w 2 Ai 
Ai w\ 



x 2 o'x), 

(po,Pi,- 
Xn \ 



(S.20) 
,p N ), and 

(S.21) 



w 



N 



J 



There is an orthogonal matrix G for which 



n 



-(p T p + x T o x), 



(S.22) 



where X = G T X, P = G T P, and is diagonal with 
&j = (0 )jj being an eigenvalue of Q 2 . We calculate 



N 

i\x (t)x \0) = GojGoj'iOlxj^Xj'lO) 

j,3'=0 



N 



r 2 



itVLj 



and obtain with that 



N 



Gl 



S{u) = 2ttw -^*( 

3=0 



OJ - 



(S.23) 



(S.24) 



(S.25) 



= 4e(w)wolm[ft(0 2 ,u/ 2 - »0 + )o,o]- (S.26) 

In the last line, the matrix element of the resolvent 
7?-(0 2 ,w 2 ) = (w 2 — H 2 ) -1 is to be taken in the basis in 
that $7 2 has the form of Eq. (S.21). It can be calculated 



following Ref. |4J. The result is 

4e(u)uo[(F h F h ) u /2 + 0+] 



S(w) 



[a; 2 -wl -2x(^ 2 )] 2 + [0++ (F h F h ) u /2}^ 

(S.27) 

MF h F h ) a 



l 

27 



dVL- 



n 2 



(S.28) 



Note that S(u>)\ g= o = 2tt5(uj — u>o). However, the spec- 
trum of any realistic microwave resonator at g = will 
be a Lorentzian with full linewidth k at half maximum. 
We use the case g = to relate the so far infinitesimal 



real number + in Eq. (S.29) to n by demanding 



49(cj)^o0+ 



(w 2 -u; 2 ) 2 + (0+) 2 (w - wo) 2 + («/2) 2 



(S.31) 



For ^ <C c^o, it is sufficient to focus on the vicinity of the 
strongly pronounced peak of S(u)\ g= o ai u = uo (i.e., 
on cj — uo <C cJo) ? and we find that here Eq. (S.31) is 
fulfilled for + = kujq. Inserting this expression in Eq. 
flS.29| ) finally leads to Eq. @ of the main text. 



Fig. [S2] complements Figs. |2|b) and [2jc) of the main 
text by showing S{uo) in the limiting cases g/J <C 1 [Fig. 
^2ja)] and #/J > 1 [Fig.[S2jb)]. In Fig.^a), we choose 
the parameters J /uoq and k/ujq as in Fig. 2(b), but g/ujQ = 
0.05. Where the Ising chain is off-resonant with cjo, the 
spectrum is qualitatively similar to the one of Fig. |2|b). 
Also here one observes the dispersive shift (oc g 2 ) of the 
resonator frequency in analogy to the N = 1 case and a 
broad side maximum of width ~ 4 J (blue and green lines 
in the inset). Both are less pronounced than in Fig.[2|b) 
due to the lower value of g. On resonance (£ « cjo/2J), 
though, the double peak structure reminiscent of the TV = 
1 case is no longer visible. Instead, S(u) is a Lorentzian 
around ujq with full width at half maximum given by 
2g 2 / J (as long as n <C g 2 / J and uo is within the band of 
the Ising chain). Indeed, assuming small g/J, one may 
replace p(u) by its maximum 2/ J and take x(^ 2 ) ~ in 
Eq. ( |S.29| ). One can then verify 



S(lj)\ 



2 9 yj 



(S.32) 



( W - Wo ) 2 + (5 2 /J) 2 ' 

In Fig. |S2[b ), we choose the parameters g/uJo and k/ujq 
as in Fig. |2[b), but J/ujq = 0.05. This case has already 
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transverse field £ = SI/2J transverse field £ = SI/2J 



FIG. S2. Spectrum S of a resonator coupled to an Ising chain 
(N — > oo) vs. frequency uj and normalized transverse field £. 
(a) The case g/ J (the parameters are g = 0.05, J — 0.1, 
and ft = 10 -4 ). Inset, 5(a;) for £ = 3.9,5,6.1 (blue, red, 
green), (b) The case g/J^>l (the parameters are g = 0.12, 
J = 0.05, and k = 10" 4 ). Inset, 5(w) for £ = 7.8,10,12.2 
(blue, red, green). All parameters are measured in units of 
ujo. The dashed lines are the first two excitation energies of 
7i for the same parameters, but N = 1. For better visibility 
of the features, values of S(uj) > 1 [S(uj) > 8] in the density 
plot of (a) [(b)] are plotted white. The lines in the insets 
correspond to cuts along the arrows in the main plots. 



much similarity with the usual single-qubit case. Off res- 
onance, the resonator experiences again the same disper- 
sive shift as for N = 1. On resonance, the broad double 
peak structure of Figs. [2|b,c) with width 4J has devel- 
oped into two sharp Lorentzians separated by ~ 2g as for 
N = 1 (red line in the inset). The chain is visible only 
as faint band of width 4J in between these peaks. 

In order to illustrate finite-size effects on the res- 
onator spectrum S(lj), we calculate the spectrum p(uS) 
of a finite transverse-field Ising chain. It is given by 
the Fourier transform of Eq. (S.14) and reads p(uj) — 
27r^ /c 0| ± 5(uj — Afc). We assume that the delta peaks 
in p are broadened by decay processes and replace them 
with Lorentzians centered around A^ and having a full 
width at half maximum of 7. Together with Eq. Q of 
the main text, this yields the spectrum S(uj) of a res- 
onator coupled to a TFIC of finite length. In Fig. |S3[ 
we plot S(uj) for similar system parameters as in Fig. [2] 
of the main text (g/J « 1), but N = 20. Signatures of 
the QPT at £ = 1, the dispersive shift of the resonator 
frequency, and the double-peak structure on resonance 
with 4 J separation of the peaks (rather than 2g as in the 
case N — 1) are present also for N = 20. We remark that 
compared to the case N — >> 00 (Fig. [2]), the ratio g/J has 
to be slightly increased for N = 20 (Fig. S3) such that 
the double peak structure of S (lj) on resonance is clearly 
visible. This is because the weight of the edges of the 
band of the Ising chain in the spectrum p(uj) increases 
with N. 
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FIG. S3. (a) Spectrum S of a resonator coupled to a fi- 
nite Ising chain (N — 20) vs. frequency uo and normalized 
transverse field £. The parameters are g — 0.12, J — 0.08, 
k — 10~ 4 , and 7 = 5 x 10~ 3 (in units of ljq). (b) S(u) for 
£ = 6.1. This curve corresponds to a cut along the arrows in 
(a). 

PROPAGATION OF A LOCALIZED 
EXCITATION IN THE ISING CHAIN 

This section contains the explicit evaluation of the 
RHS of Eq. ([5) of the main text. Further, it is shown 
that by deliberately detuning the transition frequency of 
one qubit, the effective length of the TFIC can be modi- 
fied. With 



k 



(S.33) 
(S.34) 



where <pk,j and ipj^j are determined by Eqs. (S.5) [and 
explicitly given in Eqs. (S.9) and (S.10)], we reformulate 
Eq. ([5]) in terms of fermions, 



<ai>(t) = <0|LiM J -(t)L i (t)Li|0>. 



(S.35) 



The RHS of this equation can be evaluated using Wick's 
theorem, which was first used in this context in Ref. pQ. 
One finds 



j4>k,j + £ e " v " h ' [Vk.lVk'A 



<"2>(*)=-£^ 
k 

x (^k,j4>k',j + ipk',j<t>k,j)]- 



(A fc -A fc /)tf 



k,k f 



(S.36) 



This formula was used for the plots in Fig. [3] of the main 
text. 

If the transition frequencies Qj of the qubits can be 
tuned individually, one can intentionally detune one 
qubit from the rest of the chain and observe how the 
system dynamics changes depending on the detuning. 
Fig. S4 shows the time evolution of (cr J z )(t) after a lo- 
cal excitation has been created on the first site for the 
same system parameters as in Fig. [3j but with qubit 11 
strongly detuned from the others, explicitly fin = 1.3Qj 
for j 7^ 11. This local inhomogeneity acts as a barrier 
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FIG. S4. Nonequilibrium time evolution of (a J z ) after a ty- 
pulse on the first qubit in a transverse-field Ising chain of 
length N — 20 in the paramagnetic phase (normalized trans- 
verse field £ = 8). Qubit 11 is strongly detuned from the rest 
of the chain. Values > —0.5 are plotted white. The measur- 
able observable is singled out left. 

for the propagating excitation and leads to its reflection. 
The revival of the measurable observable (cr\)(t) takes 
place at t « N/2J rather than at t w N/J as in Fig. [3] of 
the main text. Thus, strongly detuning one qubit from 
the others effectively changes the length of the chain. 



QUENCH DYNAMICS OF THE 
MAGNETIZATION AND THE END-TO-END 
CORRELATIONS 

We calculate the time evolution of (a J z ) and (cr^cr^) 
that follows a sudden change from £ = £ a to £ = £& at 
t = 0. We plot and discuss the result for (cr^cr^) and 
provide a plot of (a J z ) in addition to Fig. [Z] of the main 
text. In the following, quantities belonging to Hi, a are 
labelled by a (like AjJ), and analogously for 

First, we focus on 



= a(0\ 



e" Ul > bt (jie 



■iUi 



b£ |0) a . 



(S.37) 



To evaluate the RHS, we use the usual mapping to free 
fermions [U[2]: We express a J z by r] k and rf£ whose time 
dependence is trivial. Then we express these operators 
by r] k and rf£ whose action on |0) a is known. One obtains 



EC-C 



k,k' 



[X a k: b k ,cost(A b k + A b k ,) + cost(A> - A b k ,)}}. (S.38) 
Here, 



X?' 



k,kf 



U) T H« + (h b ) T G a ] [(G") T g b , + (H"fh b ,], 
Y k % = HfH a + (h b k ) T G«] [(H«fg b k , + (G a ) T h b ,}, 

and G and H are matrices containing the g k and hk 
as columns, respectively. Complementary to Fig. [4] of 
the main text, we plot in Fig. S5 (<j J z )(t) for a quench 
£ a = 8 — >• £5 = 1 in a transverse-field Ising chain with 




qubit number j 



FIG. S5. Time evolution of (a J z ) in a transverse- field Ising 
chain of length N = 100 after a quench of the normalized 
transverse field £ = 8 — > 1. Values < —0.55 (> —0.52) are 
plotted black (white). 

length N = 100. We focus here on t < T and choose a 
relatively large chain to strongly contrast the initial ap- 
proach of (<j 3 z ) to a constant value with the effects of the 
finite system size. The choice of the non-generic value 
£& = 1 minimizes dispersion of and, thus, of the ve- 
locities of the quasiparticles. The features of (cr J z )(t) for 
t < T described in the main text are more pronounced 
and clearly visible in Fig. [S5| 

Let us now turn to the quench dynamics of the end- 
to-end correlator 



(ala»)(t)= o (0|e 



m 



|0)„. (S.39) 



We remark that similar quantities have been studied in 



|5J. To evaluate the RHS of (S.39), we use that T~Lj com- 
mutes with e t7V ^k=i c k c k f or a n Consequently, |0) a is 
also an eigenstate of the latter operator pQ. It is now 
easy to see that 



[cUt) 



C N (t)]( 

c N {t)]\0) a , 



: |0) a 



(S.40) 
(S.41) 



where 0(t) = e mi ^Oe- iHl ^ for an operator O. The 
same strategy as for {a J z ) leads to 

<^)(t) = E<i<^ +2 E^M^'^ x 



k,k' 



[X££cosi(A b fe + At) - Y k a : b lCOS t(A b k - A b k ,)]}, (S.42) 



with X k ' k/ and Y kk/ defined above. This result is plotted 
in Fig. [S6] for a quench within the paramagnetic phase. 
The observable (cr^a^ ) is an order parameter of the Ising 
chain in equilibrium and does not develop a nonzero mean 
value for quenches within the paramagnetic phase. How- 
ever, at t « N/2vq = T/2, where vq = 2J in the para- 
magnetic phase, oscillations of (cr^a^) arise. After an 
abrupt increase, their amplitude decreases again, and this 
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FIG. S6. Time evolution of the end-to-end correlator (cr^cTx) 
in a transverse-field Ising chain of length N — 30 after a 
quench of the normalized transverse field £ = 8 —> 1.5. 

pattern quasiperiodically repeats with period T = N/vq. 
The observed behaviour of the end-to-end correlator can, 
again, be understood in the QP picture (see |5J for a re- 
lated analysis). Among the pairs of momentum-inverted 
QP trajectories with the same origin only those origi- 
nating at j = N/2 have trajectories hitting the system 
boundaries simultaneously. Since only contiguously gen- 
erated QPs carry quantum correlations, only the QPs 
generated at j = N/2 can build up correlations between 
the surface spins which will manifest themselves in a 
nonzero value of (crla^). These QPs arrive for the first 



time at the surface spins at t = N/2vq = T/2, are then 
reflected, and build up correlations between the surface 
spins each time they have travelled through the whole 
chain, that is, after multiples of T = N/vq. This ex- 
plains the two different time scales T/2 and T and the 
quasiperiodicity (for t > T/2) of the end-to-end correla- 
tor (crla^)(t). Slower QPs generated at j = N/2 will also 
arrive simultaneously but delayed at the surface spins. 
They are responsible for the slow decay of the oscilla- 
tions of the correlator for t > T/2 and will, for large t, 
eventually smear out the quasiperiodic structure. 
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